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Abstract 

We study an asymptotic behavior of the probability of formation of a ferromagnetic 
string (referred to as EFP) of length n in a ground state of the one-dimensional 
anisotropic XY model in a transversal magnetic field as n — > oo. We find that it is 
exponential everywhere in the phase diagram of the XY model except at the critical 
lines where the spectrum is gapless. One of those lines corresponds to the isotropic 
XY model where EFP decays in a Gaussian way, as was shown in Ref.[l]- The 
other lines are at the critical value of the magnetic field. There, we show that EFP 
is still exponential but acquires a non-trivial power-law prefactor with a universal 
exponent. 
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1 Introduction 



Emptiness Formation Probability (EFP) is a special correlation function which 
plays an important role in the theory of integrable models [2]. In the case of 
spin-i chains at zero temperature it is defined as 

in) = ' ■ ■ 



n 



where cTj^'^'^ are Pauli matrices, the integer j labels the sites of the Id lattice, 
and the average in (1) is performed over the ground state of the Hamiltonian 
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of the spin chain. EFP measures the probabiUty that n consecTitivc spin sites 
are all found aligned downward in the ground state. The asymptotic behavior 
of P(n) at large n is, probably, the most interesting one from the physical 
point of view, as it is believed that it measures the tendency of the system to 
a phase separation. 

There has been considerable progress in finding this asymptotic behavior in the 
case of the integrable XXZ spin chain in critical regime [2,3,1,4,5,6,7]. At finite 
temperature the asymptotic behavior of EFP is exponential [2] P{n) ~ e"'"'", 
while at zero temperature in critical regime it is Gaussian P{n) ~ e^"" . The 
qualitative argument in favor of a Gaussian decay of EFP at zero temperature 
was given in Rcf. [8] in the framework of a field theoretical approach. There, 
it was argued that an optimal fiuctuation of the field corresponding to EFP in 
the critical model will have (very roughly) the form of a round droplet in space- 
time with area ~ and will have an action proportional to n^. The asymptotic 
form of P(n) is defined by the action of this optimal configuration and is 
proportional to P{n) ~ e""" , which explains the Gaussian decay of EFP at 
zero temperature. Similarly, at finite temperature T in the limit n ^ 1/T 
the droplet becomes rectangular (one dimension n is replaced by the inverse 
temperature 1/T) and the action cost is proportional to n, giving P{n) ~ e~'^^. 
This argument is tied to the criticality of the theory. Indeed, for noncritical 
theory one should not expect the same scaling of time and space dimensions 
of an optimal configuration. It is, therefore, very interesting to see how the 
asymptotic behavior of EFP is related to the criticahty of a theory. 

In this paper we examine the relation between an asymptotic behavior of EFP 
and criticality using the example of the Anisotropic XY spin- 1/2 chain in a 
transverse magnetic field h 

-^E^I (2) 

j 

This model has been solved in [9] in the case of zero magnetic field and in [10] 
in the presence of a magnetic field. The phase diagram is depicted in Fig. 1. 
The phase diagram of the XY model has obvious symmetries 7 — > —7 and 
h —h. However, the latter one is broken by EFP itself. Therefore, we show 
only the part of the diagram corresponding to 7 > 0. The phase diagram has 
both critical and non-critical regimes. The first critical regime - the Isotropic 
(7 = 0) XY model - is labeled VLq. There is another critical line h = 1 which 
is labeled where excitations are gapless (and similarly line Q_ aX h — — 1). 
The rest of the phase diagram is non-critical with finite gap for excitations 
and is divided into domains E+, Sq, and S_ for /?. > 1, — 1 < /?. < 1, and 
h < —1 respectively. Fig. 1 includes the line 7 = 1 (F/) corresponding to an 
Ising model in transverse magnetic field and the line 7^ -|- /i^ = 1 (Fg) on 
which the wave function of the ground state is factorized into a product of 



i + 7 



2 




single spin states [11]. 
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2 The Emptiness Formation Probability in the XY Model 



It was shown in Refs. [9,10] that spin-spin correlation functions in the XY 
model can be represented as determinants of some Toephtz matrices so that 
one can find their asymptotic behavior using known theorems from the theory 
of Toeplitz matrices. This approach has been applied to EFP for the isotropic 
(7 = 0) XY model in magnetic field in Ref. [1]. In this section we (slightly) 
generahze the derivation of Ref. [1] to the case of 7 7^ and represent EFP as 
a determinant of a Toeplitz matrix. 

We follow the standard prescription [9] and reformulate the Hamiltonian (2) 
in terms of spinless fermions ijji by means of a Jordan- Wigner transformation: 



k<j 



aJ = 2V^>,-l, 



(3) 

(4) 



where, as usual, — {a^ ± i(j^)/2. In Fourier components -0^ = Y.qi'q^^'^^ we 
have: 



H = ^\2 (cosg - h) iplipq + i7 sing iplip-q - 17 sing ip-qip,^ 



(5) 



The Bogoliubov transformation 

Xq = cos^ i)q + i sin^ i)lq (6) 



with "rotation angle" -dq 

e^^" ^ —(cosq - h + ij sin q), (7) 

Sq 



brings the Hamiltonian (5) to a diagonal form J2q ^qXlXq with the quasiparticle 
spectrum 



£q 



\J (cos g — /i)^ -I- 7^ sin^ g. (8) 



The fermionic correlation functions are also easy to obtain from (5). In the 
thermodynamic limit they read [9,10] 
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F,, ^ = = / ^ ^ sin ^,e^''(^--'=) (9) 

— TT 

G*-<^.^J) = /|iif^e"<^-). (10) 



We express EFP as the fermionic correlator [1] 

P(n) = (n^,V^]). (11) 



and using Wick's theorem on the r.h.s of (11), we obtain 

P{n) = Pf(M) = ^det(M), (12) 



where Pf(M) is the Pfaffian of the 2n x 2n skew-symmetric matrix M of 
correlation functions 



Here F and G are n x n matrices with matrix elements given by Fj^ and Gjk 
from (9,10) respectively. 

We perform a unitary transformation 

M'^UMU'^(_»^, '")■ U^^(; -/). (14) 



where I is the n x n unit matrix and S„ = G + iF and Sn^ = G — iF. This 
allows us to calculate the determinant of M as 

det(M) = det(M') = det(S„) • det(S„^) = |det(S„) f . (15) 



The matrix Sn is a n x n Toeplitz matrix (its matrix elements depend only on 
the difference of raw and column indices [15]). The generating function a{q) 
of a Toeplitz matrix is defined by 

(S„),. = / ^ a{q)e''^^^-'^ (16) 
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and can be found from (9,10) as 

a{q) = 1 (l + e^^^) = ^ + ^ " ^ + ^7sin g ^^^^ 

2 2 2J{cosq — h)'^ + j^sin'^ q 

Thus, the problem of calculation of EFP 

P(n) = |det(S„)|, (18) 



is reduced (exactly) to the calculation of the determinant of the nxn Toephtz 
matrix Sn defined by the generating function (16,17). The representation (18) 
is exact and valid for any n. In our study we are interested in finding an 
asymptotic behavior at large n ^ oo. 

In the following section we apply known theorems on the asymptotic behavior 
of Toeplitz determinants to extract the corresponding asymptotes of P{n) at 
n — > cxD. The reader who is interested in results can skip the following section 
and go directly to Sec. 4 where we present and discuss our results. 



3 The results for the asymptotic behavior of EFP 



The asymptotic behavior of EFP at n — > oo is exactly related to the asymp- 
totic behavior of the determinant of the corresponding Toeplitz matrix (16,17,18) 
and can be extracted from known theorems and conjectures in the theory of 
Toeplitz matrices. These types of calculations have been done first in [9,10] 
for spin-spin correlation functions. It is well known that the asymptotic be- 
havior of the determinant of a Toeplitz matrix as the size of the matrix goes 
to infinity strongly depends upon the zeros of the generating function of the 
matrix and upon its phase. 

In the case when the generating function a{q) has only isolated singularities 
(zeros and phase discontinuities) the most general known result for the asymp- 
totic behavior of a Toeplitz determinant is called generalized Fisher- Hartwig 
conjecture (gFH) and is due to Basor and Tracy [15]. The result can be written 

as 

P(n) ~ ^£;in-^'e-^*", (19) 



where the decay rates jSi, the exponents Aj, and the coefficients Ei are known 
functionals of a{q) [15]. The original Fisher- Hartwig conjecture (FH) [13] con- 
sists of just one term of the sum in (19). We apply gFH to the generating 
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function (17) and extract the asymptotic behavior of P{n). All necessary the- 
orems on the asymptotic behavior of Toeplitz determinants can be found in 
Ref. [12, 13, 14, 15, 17, 16]. For the most recent review of gFH see also Ref.[18]. 

In some parts of the phase diagram we have to rely on "conjectures", i.e., 
results which are not proven yet. Therefore, it is prudent to look for some 
independent ways of checking those results. We did some numeric calculations 
on these Toeplitz determinants and the results of these calculations support 
the analytic results presented below [19]. Also, we used some hmiting cases 
where EFP or an asymptotic behavior of Toeplitz determinant can be obtained 
in an independent way to confirm our conclusions. E.g., on the line EFP 
can be calculated exactly and the asymptotic behavior of the result is the same 
as the one obtained using FH for the whole region Sq of the phase diagram. 

In the rest of this section we present the results for the asymptotic behavior of 
P(n) in different regions of the phase diagram of the XY model (Fig. 1). In this 
letter we concentrate on the calculation of A and /3(/i, 7) leaving our findings 
for the coefficient £^(/t, 7) as well as the results of our numeric calculations for 
a future, more extended publication [19]. 



3.1 The non-critical theory (Ti± and T^q) 



3.1.1 The region T,^ 

In this region (7 7^ 0,/t < —1) the generating function (17) is nonzero for all q 
and the FH conjecture is reduced to a (rigorously proven) Strong Szego Limit 
Theorem. It gives 

P(n) ~ E(/i,7)e-^('^'^)" (20) 



with 



(5{Kl) = - j'^\og\a{q)\ 




cos q — h 



cos q — h) +7^ sin^ q, 



[21] 



The integral in (21) is convergent for all h and all 7 ^ and f3{h,'~f) is a 
continuous function of its parameters. In fact, it turns out (see below) that this 
expression for (3{h,j) describes correctly the decay rate of EFP everywhere 
in the phase diagram except for 7 = 0. However, both pre-exponential factors 



7 
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Fig. 2. Plot of the decay factor /? as a function of the parameters 7 and /i. The 
function diverges for 7 = and is continuous for ^ = ±1 (although has a weak 
singularity at /i = ±1). 

and the mathematical status of the results are quite different in different parts 
of the phase diagram (see below) . 

In Figure 2, /9(/i, 7) is plotted as a function of h at several values of 7. One 
can see that /3(^, 7) is continuous but has a weak (logarithmic) singularity at 
h — ±1. This is obviously one of the effects of the criticahty of the model on 
the asymptotic behavior of EFP. 

3.1.2 The region J^o 

In Eq (7 7^ 0,-1 < h < 1) the generating function a{q) vanishes at g = ±n 
and acquires an overall phase of tt as g goes from — tt to n. The Szego theorem 
is not applicable anymore. One has to rely on the FH conjecture (which, in 
fact is proven for this case - the case of one singularity [17]) to derive the 
result, which is identical to (20,21). The expression for the coefficient £'(/i, 7) 
is different from the one in E_ but its form is beyond the scope of this paper. 

3.1.3 The region J:+ 

In (7 7^ 0, /i > 1) the generating function a{q) vanishes at the points 
q — and q = ±7r, and its phase has jumps at these points. In this case an 
application of FH leads to some ambiguity, e.g., one obtains two values for 
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P(h, j): Pi — (3 and = /3 + vri with f3 from (21). This ambiguity is resolved 
in the use of the gFH which gives EFP as a sum of two terms so that both 
values of /5 are used. 

We obtain 

P(n) - E{h, 7) [1 + 7)] e-"'^('*'^). (22) 



Thus, in this region we again have an exponential asymptotic decay of EFP 
but, the constant prefactor is different for even and odd n. 



3.1.4 The line r2.' an exact calculation 

Before proceeding to the analysis of EFP on the critical lines let us check our 
results (20,21) on the special line^ in the phase diagram defined by 

h^ + -f^^ 1. (23) 



It was shown in Ref.[ll] that on this line the ground state is a product of 
single spin states and is given by 

3 3 



cos 




T,j) + (-l)^sin 




(24) 



where | t, j) is an up-spin state at the lattice site j etc. One can directly check 
that the state (24) is an eigenstate of (2) if the value of the parameter 9 is 



cos'^^^-^ (25) 
1 + 7 



and (23) is satisfied. It is also easy to show [11] that this state is, in fact, the 
ground state of (2). 

The probability of formation of a ferromagnetic string in the state (24) is 
obviously 



We are grateful to Fabian Essler who suggested us to check our results on this 
special line and pointed out the reference [11] to us. 
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which is an exact result on the hne (23). The value of /3(/i,7) which immedi- 
ately follows from this exact result is 

„(A=;rr7,,)._iog(l-iyiHD (27) 

Indeed, one can easily check that the integral in the general formula for the 

decay rate (21) under the condition (23) gives precisely (27). In fact, the 
Toephtz matrix generated by (17) becomes triangular on the hnc (23) with 
diagonal matrix element {Sn)jj = sin^(^/2) and the determinant of Sn is 
exactly (26). 

3.2 The critical line fig (j — 0) 

The case 7 = 0, corresponding to the Isotropic XY Model, has been extensively 
studied in Ref. [1]. For 7 = the generating functions (17) is reduced to the 
one found in [1]. 

For h < 1, the generating function a{q) has a hmited support: to find the 
asymptotic behavior of the determinant of the Toeplitz matrix one can apply 
Widom's Theorem [16] and obtain [1]: 

P{n) ~ 2Ae3^'(-i)(l - /i)-in-3 ( ^ ) ' (28) 



and we see that in this case, EFP decay as a Gaussian with an additional 
power-law pre-f actor. ^ 

For \h\ > 1, the theory in not critical anymore and the ground state is com- 
pletely polarized in z direction, giving a trivial EFP P(n) = for /i > 1 and 
P{n) = 1 for /i < -1. 

3.3 The critical lines Q±(h = ±1) 

3. 3. 1 The critical line n+(h^l) 

For h = 1 the generating function a{q) vanishes ai q = ±71 and its phase has 
jumps at q — 0, n. The gFH in this case produces many terms of the form (19). 

^ Notice, that the formula (28) was in fact obtained in [20] in the context of unitary 
random matrices (which are known to be equivalent to free fermions). 
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However, in contrast to the S+ region, all terms are suppressed by power law 
factors of n with respect to the leading one. The leading term has fi from (21) 
and critical exponent A = ^: 

P{n) ~ E{^)n-^ Q-nf3{-y,i)_ (29) 

The next subleading term (faster decay at n — > oo) is given by 

£;'(7)(-l)"n-^e-"'^(i'^^ (30) 



Although inclusion of the latter (subleading) term is somewhat beyond even 
the generalized FH, we write the sum of these two terms as a conjecture for 
EFP at /i = 1 and obtain 



P(n) ~ £;(7) n~T6 i + (_i)M(7)/n2 



-n^(l,7) 



(31) 



As this result relies on this unproven conjecture, we present our numeric results 
for this case in Figure 3. Indeed, we see that the form (31) is in good agreement 
with numerics. Wc conclude that h — 1 EFP decays exponentially at n ^ 
oo but with an additional power law pre-factor and a damped oscillatory 
component. 

Remark. It is curious to notice that 1/16 and 9/16 remind of the scaling 
dimensions of spins and It looks as if EFP operator (1) among other 
things has inserted square roots of spins at the ends of the string. 



3.3.2 The critical line Q-(h^ -1) 

For h — —1 the generating function a{q) does not vanish but has a phase 
discontinuity at g = ±7r. We apply FH to obtain 

P{n) - E(7) n-Te e-"/3(-i,7). (32) 



and stretching gFH the same way as we did for /i = +1 we include a subleading 
term 

P{n) ~ ^(7) n-Te [1 + A{^)/n^ e^^^^-^'^l (33) 



This form is also in good agreement with our numeric results. 

^ See Ref.[10] where it was shown that power laws are 1/4 and 9/4 for and 
correlators respectively. 
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0.86 - 
0.84 - 
0.82 - 
0.8 - 
0.78 - 
0.76 - 




Fig. 3. Results of the numeric calculation of the Toeplitz determinant is shown as 
P(n)e^"n^/^6 vs n at 7 = 1, /i = 1. The value of /? = log 2+2G/7r with Catalan's con- 
stant G is obtained from (21). The solid line is the analytic result E{l + {—l)'^A/n2) 
with E = 0.8065..., A = 0.2399... obtained by fitting at large n. To make the plot 
more readable we show only every 11th point (for n = 1, 12, 23, . . .). 

4 Discussion and Conclusions 

An asymptotic behavior of Emptiness Formation Probability P{n) as n ^ cxd 
for the Anisotropic XY model in a transverse magnetic field as a function of 
the anisotropy 7 and the magnetic field h has been studied. We summarize 
our results in Table 1. All these results are new except for the isotropic case 
^0 (7 = 0) which was studied in Ref.[l]. 

Our main motivation was to study the relation between the criticality of the 
theory and the asymptotics of EFP. Let us look at the results on the critical 
lines and Q±. The Gaussian behavior on is in accord with the qualitative 
argument of Ref. [8]. However, on Q± the decay of EFP is exponential instead of 
Gaussian, and apparently contradicts the qualitative picture of Ref. [8] (which 
is reproduced in the introduction to this paper for the sake of completeness). 
The reason for this disagreement is that &t h — ±1 the model is critical for 
the quasipartides x defined in (6). However, we are interested in EFP for the 
"original" Jordan- Wigner fermions ip, and this correlator has clearly a different 
expression in terms of x than the simple one (11) which holds at 7 = 0. From 
the technical point of view, the difference is that in the qualitative argument 
in favor of a Gaussian decay of EFP for critical systems there is an implicit 
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Table 1 



Asymptotic behavior of EFP in the different regimes. The exponential decay rate /3 
is given by Eq. (21) for all regimes. The coefficients E, A are functions of h and 7 
as appropriate. 



EFP for an Anisotropic XY model 


Region 


7, h 


criticality 


P{n) 


Eq.# 


Theorem 




h<-l 


No 




20 


Szego 


So 


-l<h<l 


No 




20 


FH 


s+ 


h > 1 


No 


E[l + {-lYA] e-"'' 


22 


gFH 


To 




No 




26 


Exact 




h = l 


Yes 


E n-Vi6 [1 + (-l)n^/^ ] Q-nQ 


31 


gFH 


n_ 


h = -l 


Yes 


En-^/^^ [1 + A/^] e-""^ 


33 


gFH 


no 


7 = 0, -l<h<l 


Yes 




28 


Widom 



assumption that the density of fermions (or magnetization) is related in a local 
way to the field responsible for critical degrees of freedom (free boson field 0). 
This assumption is not valid on the lines h = ±1. Although the theory is 
critical on those hnes and can be described by some free field 0, the relation 
between the magnetization and this field is highly nonlocal and one can not 
apply the simple argument of [8] to the XY model at /i = ±1. 

Although EFP at the critical magnetic field does not show a Gaussian be- 
havior, there is an important difference between the asymptotic behavior of 
EFP on and off the critical lines. Namely, a power-law pre-factor appears on 
all critical lines. For the XY model it is universal (constant on a given critical 
line) and takes values A = 1/4 for 7 = [1] and A = 1/16 on the lines h = ±1. 
It would be interesting to understand which operators determine these par- 
ticular values of "scahng dimensions" of EFP (see the remark at the end of 
section 3.3.1). 

In some parts of the phase diagram we used the so-called generalized Fisher- 
Hartwig conjecture [15] which is not proven yet. However, our numeric calcu- 
lations support our analytical results, see, e.g.. Figure 3. We note that to the 
best of our knowledge this is the first physically motivated example where the 
original Fisher-Hartwig conjecture fails and its extended version is necessary. 
We also suggested that the gFH could be used to find the subleading correc- 
tions to the asymptotic behavior, as we did for h — ±1 in (31,33). Prom the 
physical point of view, the oscillations of EFP that we obtain for /i > 1 from 
the gFH are related to a spontaneous Z2 symmetry breaking which occurs at 

In conclusion, we notice that it is straightforward to generalize our results 
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for nonzero tcmpcratTirc. The only modification is that at T 7^ the thermal 
correlation functions must be used instead of (9,10). Then, the generating 
function (17) is non-singular everywhere and we have an exponential decay of 
P{n) in the whole phase diagram according to the standard Szego Theorem. 
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